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ABSTRACT 


How much electromagnetic energy crosses the photosphere in evolving solar 
active regions? With the advent of high-cadence vector magnetic held observa¬ 
tions, addressing this fundamental question has become tractable. In this paper, 
we apply the “PTD-Doppler-FLCT-Ideal” (PDFI) electric held inversion tech¬ 


nique of Kazachenko et al. (2014) to a 6-day HMI/SDO vector magnetogram 


and Doppler velocity sequence, to hnd the electric held and Poynting hux evolu¬ 
tion in active region NOAA 11158, which produced an X2.2 hare early on 2011 
February 15. We hnd photospheric electric helds ranging up to 2 V/cm. The 
Poynting huxes range from [—0.6 to 2.3] x 10^° ergs-cm“^s“^, mostly positive, 
with the largest contribution to the energy budget in the range of [10®-10^°] 
ergs-cm“^s“^. Integrating the instantaneous energy hux over space and time, we 
hnd that the total magnetic energy accumulated above the photosphere from the 
initial emergence to the moment before the X2.2 hare to be i? = 10.6 x 10^^ ergs, 
which is partitioned as 2.0 and 8.6 x 10^^ ergs, respectively, between free and 
potential energies. Those estimates are consistent with estimates from prehare 
non-linear force-free held (NLFFF) extrapolations and the Minimum Current 
Corona estimates (MCC), in spite of our very diherent approach. This study 
of photospheric electric helds demonstrates the potential of the PDFI approach 
for estimating Poynting huxes and opens the door to more quantitative studies 
of the solar photosphere and more realistic data-driven simulations of coronal 
magnetic held evolution. 


Subject headings: Sun: magnetic field, Sun: flares. Sun: sunspots 
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1. INTRODUCTION 

The advent of high-cadence, large-scale vector magnetic field and Doppler velocity 
measurements from instruments such as the Helioseismic and Magnetic Imager (HMI, 


Schou et ah (2012)) on NASA’s Solar Dynamics Observatory (SDO) satellite (Pesnell et ah 


2012), the Spectropolarimeter instrument (SP; Lites et ah 2013) on the Solar Optical 


Telescope (Tsuneta et ah 2008) aboard the Hinode satellite (Kosugi et ah 2007), and 


improved capabilities of ground-based instruments, such as SOLIS (e.g., Keller et ah 2003), 


make the estimation of electric helds in the solar photosphere possible. The calculation 
of the electric held from magnetic and Doppler data is critically important for various 
quantitative studies of the solar atmosphere. First, if we know both electric and magnetic 
held vectors in the photosphere, we can estimate both the Poynting hux of magnetic energy 
and the hux of relative magnetic helicity entering the corona. Second, as demonstrated 


in a magneto-frictional model by Cheung and DeRosa (2012), the ability to compute the 


electric held enables the driving of time-dependent simulations of the coronal magnetic 
held from photospheric magnetogram sequences. Combining electric held estimation with 
a magneto-frictional model of the evolving solar corona is the goal of the Coronal Global 
Evolutionary Model (CGEM) project (Fisher et al. 2015)]^ 


Kazachenko et al. (2014) modihed and extended the electric held inversion methods 


introduced by Fisher et al. (2010, 2012), to create a comprehensive technique for 


calculating photospheric electric helds from vector magnetogram sequences. The new 
method, which we dubbed the PDFI (an abbreviation for Poloidal-Toroidal Decomposition 
[PTD]-Doppler-Fourier Local Gorrelation Tracking [FLGTj-Ideal) technique, incorporates 
Doppler velocities from non-normal viewing angles (which are relevant to most solar 
observations) and a faster and a more robust Poisson equation solver, for obtaining the 


^http://cgem.stanford.edu 
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PTD solutions. After systematic, quantitative tests of the accuracy and robustness of the 


PDFI technique, using synthetic data from anelastic MHD (ANMHD) simulations (Abbett 


et al. 2004), we found that the PDFI method has less than 1% error in the total Poynting 


flux and a 10% error in the helicity flux rate if we reconstruct it at the normal viewing angle 
(6* = 0) and less than 25% and 10% errors respectively at large viewing angles (6* = 60°) 


(Kazachenko et al. 2014). 


In this paper, we take the next step forward, and apply the PDFI technique to 
observations. The flare-productive active region (AR) NOAA 11158 was observed by HMI 
nearly continuously for a six-day period over 10-16 February 2011, starting from its initial 
emergence near 14:00 UT on 10 February. We use the sequence of magnetic and Doppler 
held measurements of NOAA 11158 to derive the temporal evolution of electric held, 
Poynting, and helicity huxes during these six days. The evolution of this AR included two 
large bipoles emerging in close proximity, with strong shearing motion between the central 
sunspots dSchrijver et al. 2011 Sun et al. 2012). Over a six-day period, the AR hosted an 
X2.2 hare (with the GOES soft X-ray hux peaking at 01:56 UT on February 15) leading 
to a pronounced halo CME, three M-class hares, and over twenty C-class hares. Since this 
active region was the hrst one for which the HMI vector magnetic held data were widely 
distributed to the scientihc community, its magnetic held has been thoroughly studied: 


the fast sunspot rotation from 20 hours before to 1 hour after the X2.2 hare (Jiang et al. 


2012 Vemareddy et al. 2012a), the hare related enhancement in the horizontal magnetic 


held along the magnetic polarity inversion line (PIL; Gosain (2012); Wang et al. (2012); 


Liu et al. (2012a)), abrupt changes in the vertical Lorentz force vectors (Petrie 2012 


Alvarado-Gomez et al. 

2012 

) and horizontal Lorentz forces ( 

Petrie 

2013 

Wang et al. 

2014), 


the evolution of relative and current magnetic helicities (Jing et al. 2012), the injection of 


oppositely signed helicity through the photosphere (Dalmasse et al. 2013), the subsurface 


three-dimensional magnetic structure (Ghintzoglou and Zhang 


2013), the magnetic and 
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velocity field transients driven by the flare (Manrya et al. 2012). Nnmerons approaches 


have been nsed to calculate the energy associated with X2.2 flare and the AR as a whole: 


the DAVE4VM method (Liu and Schuck 2012; Tziotziou et al. 2013), non-linear force 


free extrapolation (Sun et al. 2012 Tziotziou et al. 2013), the Minimum Current Corona 


Model (Tarr et al. 2013), and the coronal forward-htting method (Aschwanden et al. 2014 


Malanushenko et al. 2014). In this paper, we apply the PDFI technique (Kazachenko et al. 


2014) to derive electric helds at the photosphere, and use those to estimate the photospheric 


energy and helicities fluxes’ evolution in NOAA 11158. 


The paper is structured as follows. In Section]^ we briefly review the PDFI method itself 
and its recent improvements. In Section we describe the observations of the emerging, 
flaring NOAA 11158, and quantify how the uncertainties in the HMI observations propagate 
into derived electric held and Poynting huxes. In Section we describe derived electric 
helds, Poynting, and helicity huxes in NOAA 11158, and in Section we discuss results 
and draw conclusions. In addition, in Appendix]^ we show how our Mercator-reprojected 
magnetic helds in Cartesian coordinates can be scaled to apply to relatively small regions 
on the surface of the Sun to derive the electric helds in a local Cartesian coordinate system. 
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2. METHODOLOGY: PDFI TECHNIQUE, POYNTING & HELICITY 

FLUXES 


The PDFI technique uses the evolution of the vector magnetic held B and horizontal 
and Doppler velocities V to estimate electric helds in the solar photosphere. It is described 


in detail in §2 of Kazachenko et ah (2014). The PDFI method combines the inductive 
contribution to the electric held E from solution to Faraday’s law using the Poloidal- 


Toroidal-Decomposition (PTD) technique (Fisher et al. 2010, 2012), with non-inductive 
contributions from (—V x B). To hud the velocity V, we use Doppler measurements and 
the Fourier Local Correlation Tracking (FLCT) technique. The PDFI technique is tested and 


its accuracy is characterized in detail in §4 of Kazachenko et al. (2014). In this section we 
briehy describe the basics of the PDFI technique. 


The fundamental idea of the PTD part of the PDFI method is that the magnetic held, 
B, dehned on the photospheric surface, has a solenoidal nature and thus can be specihed 
by two scalar functions B and J': 


B = VxVxi3z + Vx Jz, 


( 1 ) 


where 5 points upward from the photosphere. Taking a partial time derivative of 
Equation (j^, and demanding that B obeys Faraday’s law, 

— = -(V X cE), (2) 

we hnd a solution for the inductive part of the electric held E^ (where “P” stands for 
PTD), in terms of the partial time derivatives of the poloidal and toroidal potentials, B and 
respectively: 

cE^ = -V X i3z - jz. (3) 


As described in Kazachenko et al. (2014), solving two-dimensional Poisson equations in the 


domain, where we observe the vector magnetic helds, we determine B, J and Note 
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that the vector magnetic held data completely specify the source terms of these Poisson 


equations (see §2.1 of Kazachenko et ah (2014)). 


The total electric held E is a combination of the inductive and non-inductive parts, 


cE = —V X Bz — jz — = cE^ 


VV'. 


( 4 ) 


inductive non-inductive 

The non-inductive components to the scalar-potential part of the solution include three 
separate contributions: (1) —V x B from Doppler measurements (the “D” in PDFI), (2) 
-V X B from Fourier Local Correlation Tracking results (the “F” in PDFI), and (3) a 
scalar potential contribution added to impose the constraint E ■ B = 0, consistent with the 
ideal MHD Ohm’s law (the “I” in PDFI). When adding the —V x B contributions above, 
any inductive contributions from these terms are removed, since all inductive contributions 
are already included in the E^ solution. Our approach for handling these non-inductive 


contributions is described in detail in Kazachenko et ah (2014). 


To calculate the hux of electromagnetic energy at the photosphere, given by the 
Poynting flux vector 

S = ^(E X B), (5) 

we use the observed magnetic field vector and the electric field vector derived using the 
PDFI method. Since we are interested in the amount of energy flowing into and out of the 
corona, we focus most of our attention on the vertical component of Poynting flux. 


5*2 = ^ {E^By - EyB^) . ( 6 ) 

This depends upon the horizontal components of both the electric field and the magnetic 
held. We further decompose Sz into two contributions, the hux of potential-held energy, 
and the hux of free magnetic energy. The basic idea is that the horizontal magnetic 
held B/i can be divided into a potential-held contribution B^, and a contribution, B{, 








due to currents that flow into the atmosphere from the photosphere (Welsch 2006| : 

X ( . In this paper we use the Green’s function to And the potential 

held contribution, and subtract this contribution from the measured horizontal fields to And 
B{. More discussion of the potential and free energy decomposition can be found in §3.1 of 


Kazachenko et ah (2014). 


To calculate helicity flux rates, we use Equation (62) from Berger and Field (1984): 
'dHn 


dt 


-2 j (A^ X E) ■ z da = -2 j {A^Ey - A^E^) da, 


( 7 ) 


PDFI 


where ) = V x B^z is the vector potential that generates the potential 


held B^ in volume V above the photosphere, which matches the photospheric normal held 


Bz aX. z = 0. (Note that a similar expression for helicity. Equation (41) in [Kazachenko 
et ah (2014), contains a typographical error - Ez should be replaced by E^.) Here B^ can 


be found by solving the Poisson equation A5 in Fisher et ah (2010|. Adopting the naming 


convention from §2.3.4 in Kazachenko et ah (2014), the total helicity flux rate derived from 


'PDFI 


the PDFI electric held E will be referred to as the PDFI helicity flux rate or 

If we have an ideal electric held cE = —V x B, then the helicity flux rate becomes 


(Berger 1984): 


dH 


R 


dt 


= -2j [(A^ ■ Vh)H, - (A^ • B^)Vz\da. 


( 8 ) 


DFI 


For observations near disk center, where the line-of-sight direction approximates the vertical 
direction, in Equation ([^, we can make the assumption that Vh can be determined with 
our FLCT flow estimates, and 14 from our Doppler velocity measurements. Adopting the 


naming convention from §2.3.4 in Kazachenko et ah (2014), the total helicity derived this 


way would correspond to DFI electric held solution (Doppler FLCT Ideal). 


When comparing helicity fluxes calculated using the PDFI and DFI techniques, it is 
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important to remember that PDFI and DFI methods are independent of each other, hence 
their results are not necessarily consistent. For more details on PDFI and DFI helicity fluxes 


and the quality of their reconstruction using ANMHD simulations, see §3.2 in Kazachenko 


et al. (2014). 
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3. DATA REDUCTION: NOAA 11158 


We derive the evolution of the magnetic held, electric held, Poynting, and helicity 
huxes in NOAA 11158 using series of HMI vector magnetograms and the PDFI method. 

A six-day, uninterrupted, 12-minute-cadence data set allowed us to study in detail both 
long-term, gradual evolution, as well as rapid changes centered around the region’s X-class 
hare. In this section we describe the data set and the coordinate system re-projection we 
use. 


NOAA 11158 was the source of an X2.2 hare on 2011/02/15, starting at 01:44 UT, 
peaking at 01:56 UT, and ending at 02:06 UT. A front-side halo CME accompanied the 


hare (Schrijver et ah 2011). Prior to the X2.2 hare, the largest hare in this region was an 


M6.6 on 2011/02/13 at 17:28 UT, a little more than 30 hours before the X2.2 hare. 

HMI observed AR 11158 in great detail, routinely generating hltergrams in six 
polarization states at six wavelengths on the Fe I 617.3 nm spectral line. From these 
hltergrams, images for the Stokes parameters, I, Q, U, and V were derived which, using 


the Very Fast Inversion of the Stokes Algorithm (VFISV) code (Borrero et al. 2011), were 


inverted into the magnetic held vector components. To resolve the 180° azimuthal held 


ambiguity the “minimum energy” method (Metcalf et al. 1994 Leka et al. 2009) was used. 


In addition, we hipped the azimuths of the transverse magnetic held vectors in all pixels 
which exhibited single-frame huctuations in azimuth of larger than 120° and for which such 


hipping reduced time variation in the azimuth (Welsch et al. 2013) 


3.1. Deriving Magnetic Vector Fields: B 

To study prehare photospheric magnetic evolution, and to baseline this evolution 
against posthare evolution, we obtained 153 hours of 12-minute-cadence 0.5"-pixel HMI 
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vector-magnetogram data, from the beginning of the active region emergence, about four 
days before the X2.2 flare, to two days after the flare: tgtart = February 10 2011 14:00 UT, 
with the active region centered at (S19, E50); and tend = February 16 2011 23:48 UT, with 
the active region at (S21, W37). 


We rotated the active region to disk center and transformed it to a local Mercator 


re-projected Cartesian coordinate system (Welsch et al. 2013). To do so, in the hrst 


step, we re-projected the observed magnetic vectors’ components onto radial/horizontal 
coordinate axes. We then converted the Cartesian output grid’s points into plane-of-sky 
(PCS) coordinates, and interpolated the radial and horizontal components of the magnetic 
held, Br and Bh, onto the remapped output grid points. Finally, to account for small, 
whole-frame shifts of the AR’s structure between successive measurements, we co-aligned 


the data to sub-pixel scale. Note that since FLCT (Fisher and Welsch 2008), and indeed 


any method of estimating the optical how (e.g., Schuck (2006)), depends upon image 


structure (e.g., gradients), conformal mappings are preferred since they are shape-preserving 


for inhnitesimally small objects. Accordingly, we use Mercator re-projection (Welsch 


et al. 2009), with equally spaced grid points as an input for FLCT. After re-projection, to 


preserve physical quantities of magnetic helds and velocities, we corrected the hux densities 
for the distortion of pixel scale introduced by re-projection; the details of the applied 
correction-factors are given in Appendix [A} 

For the minimum magnetic held to consider in the PTD, we chose a threshold of 
|B| = 250 Gauss, consistent with the upper limit of the uncertainty in the horizontal and 


vertical components of the magnetic held (Hoeksema et al. 2014). To avoid spurious signals 


in electric helds, we apply a mask, where we set any pixel’s magnetic held components to 
zero if in any of three consecutive frames it has |B| < 250 Gauss. To increase the accuracy 
of the calculated electric helds we also added a boundary area of 55-pixels width/height 


















- 10 - 


padded with zeroes around the periphery of the magnetogram (see §2.1 of Kazachenko 


et ah (2014)). The hnal data cube, after re-projection and boundary padding, consists of 


770 time steps {dt = 720 sec) and has a held of view of 665 x 645 pixels with a pixel size 
of 360.16 km, which at disk center is equivalent to the original 0.5” size of HMI’s pixels. 
Further detail on the data cube preparation and calibration, for a shorter time period, can 


be found in Welsch et ah (2013) 


Figure shows the hnal vertical magnetic held in a subregion of the full-disk data array 
after re-projection in the beginning (Panel A), middle {B,C) and the posthare {D) times 
of the magnetogram sequence. Note that the positive and negative vertical magnetic huxes, 
Bz{t), shown on the right panel, nearly balance each other; the signed magnetic huxes grow 
from essentially zero to roughly 1.4 x 10^^ Mx ( —1.4 x 10^^ Mx) at the time of the hare 
(vertical dashed line). The hux emerges in two phases - an initial, gradual phase is followed 


by a much more rapid phase, a pattern seen in the emergence of many active regions (^Fu 


and Welsch 2015). 


3.2. Deriving Velocity Vector Fields: V 


To derive the three-component velocity vector of the magnetized plasma, we used 


the following two methods: the method of Welsch et ah (2013) for calibrated line-of-sight 


(LOS) Doppler velocity component Vlqs, and FLCT for the horizontal velocities I 4 . To 
derive and calibrate the Doppler velocity for instrumental ehects and convective blueshift, 
we used three successive vector magnetograms and one Dopplergram coincident with the 


central magnetogram. Following the idea of Welsch et al. (2013), that the Doppler shifts 


measured along polarity inversion lines (PILs) of the LOS magnetic held determine one 
component of the velocity perpendicular to the magnetic held, we calibrated the quantity 
Vlos by subtracting the median Doppler velocity among all pixels on LOS PILs. To hnd 
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Fig. 1.— Panels A-D: HMI vertical magnetic field {B^) maps at 4 different times of NOAA 
11158 evolution. Panel E\ evolution of the positive and negative vertical magnetic fluxes 
of the 6-day interval, with the diamonds indicating the times of the four images (A-D) on 
the left. An X-class flare occurred at the time corresponding to the vertical dashed line. 
The black box in Panel D indicates the field of view of Figures |4]-[^ however all energy and 
helicity estimates described further are for the entire field of view. 

the horizontal velocity 14 , we determined local displacements of magnetic flux between two 
successive images in the neighborhood of each pixel, employing the following three steps. 
First, we masked the initial and hnal images with a Gaussian windowing function with an 
e-folding width of apicT = 5 pixels; second, we cross-correlated the two masked images; 
hnally, we found the peak of the cross correlation function. The vector displacement of this 
peak from zero is the inferred spatial displacement of the pattern in the neighborhood of 
the windowing function’s center. 

To calculate electric fields, and hence Poynting and helicity fluxes (Equations (|^ 
and 0 ). we use Doppler and FLCT velocities as an input into the PDF I inversion. For 
comparison, apart from the PDFI, we also use Doppler and FLCT velocities on their own as 
an independent estimate for the helicity flux rate (see Equation (|^). 

To summarize, as a result of the data reduction we obtained a six-day data cube 
consisting of 768 frames (two frames less than the original dataset to obtain FLCT velocities). 
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each of which contains data for three components of the velocity field and three components 
of the magnetic field with a field of view of 665 x 645 pixels, a pixel size of 360.16 km, and 
a time step of dt = 720 sec. 
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3.3. How do Errors in Magnetic Field Measurements Affect our Electric 

Fields Estimates? 


When testing the PDFI electric held inversion technique using MHD simulations 


(Kazachenko et ah 2014), we have the advantage of knowing that the errors in the input 


data are zero, but such is not the case with the HMI data (Liu et ah 2012b). Hence to 
estimate the uncertainties in the derived electric held and Poynting huxes, we have to 
account for uncertainties in the HMI input data, which arise primarily from estimation and 
inversion of the Stokes prohles (Liu et ah 2012b| Hoeksema et ah 2014| ). 

Fitting the core of magnetic-held-values distribution in the weak-held regions with a 
Gaussian and assuming that its width indicates the noise level in the measurements, we 
estimate the errors in Bx, By and B^ in AR 11158 dataset during six days of AR evolution 
(see Figure |^. The huctuation of the error in B^, By and B^ varies within 100 Gauss for 
the horizontal magnetic held and within 30 Gauss for vertical magnetic held due largely 


to ehects arising from SDO’s 24-hour orbital period (Hoeksema et ah 2014). We use those 
values, i.e., [100,100,30] Gauss, respectively, as noise thresholds for magnetic held values. 



Fig. 2.— Estimated noise levels in B^, By and B^, as functions of time in February 2011. 
Note the periodic variation, believed to be a function of orbital phase of the SDO satellite. 
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Fig. 3.— Pixel-by-pixel scatter plots showing the ensemble of pertnrbed electric held soln- 
tions as a fnnction of the nnperturbed electric held solntion, for all three components of E, 
for one particular realization of pseudo-random noise. 

We use Monte Carlo simulations to estimate the errors in the PDFI electric helds caused 
by the uncertainty in the HMI data. The sensitivity of PDFI electric held solutions to vector 
magnetogram noise is exacerbated by the fact that source terms of the Poisson equations 
that are solved as part of the electric held inversions involve temporal and/or spatial 
derivatives, which greatly amplify the noise. This is, however, ameliorated by the fact that 
solutions to Poisson’s equation tend to smooth out the ehect of noisy source terms. Since 
the entire inversion procedure is quite complex, our approach is to start from a given set 
of input magnetic held data, and then add pseudo-random, Gaussian perturbations to the 
data, consistent with the noise thresholds given above ([100,100,30] Gauss). By comparing 
the ensemble of electric held solutions that are perturbed about our initial solution, we can 
characterize the resulting errors of the electric held inversion. Figure shows three scatter 
plots between the original and the perturbed the electric helds for one particular realization 
of pseudo-random noise, one plot for each component. 

The resulting uncertainties in Ex and Ey are smaller than those for rehecting the 
fact that the inductive horizontal electric held from PTD, depends on 5^, which has 
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much smaller intrinsic errors than and By^ while the inductive part of the vertical electric 
field, , depends on the more-noisy B^ and By. Thus our study not only shows that the 
electric field inversion errors are not too large, but the distribution of these errors about the 
diagonal line provides quantitative estimates for the random errors, and deviations from 
unity slope provide some information about more systematic errors. To scale the deviations 
relative to the unperturbed electric fields, we fit the difference between the electric field 
from the perturbed magnetic field and the electric field from the unperturbed magnetic 
field with a Gaussian and then normalize it by dividing it by the standard deviation of the 
unperturbed data. We find the standard deviation a = [0.14,0.13,0.18,0.14] for Ey, E^ 
and Sz, respectively. 


To summarize, the levels of the errors in the HMI data lead to errors in the electric field 
up to 14% and 13% in E^ and Ey and 18% in E^. The error in the vertical Poynting flux, 
Sz, is 14%. If we add these uncertainties from the data (in quadrature) to the errors in the 
total Poynting flux that we found when testing the PDFI method (where the discrepancy, 
5Sz, ranged from 1% to 25% for viewing angles ranging from 0 to 60°, see Figure 10 of 


Kazachenko et al. (2014)), then we end up with the error ranging from 14% to 29% in the 


vertical Poynting flux. 
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4. RESULTS 


To describe the photospheric electric fields E and energy finxes S in evolving NOAA 


11158, we nse two approaches. In Section 4.1 we show the spatial distribntion of B, V, E 
and S at two times, before and after the X2.2 flare. In Section 4.2| we analyze their temporal 
spatially integrated evolntion over the six days of observations. These approaches allow ns 
to captnre both the spatial strnctnre at qniet times before the flare with the changes dnring 
the flare, and the overall long-term behavior of the active region. 


4.1. Properties Of the Active Region NOAA 11158 Before and After the X2.2 

Flare 


As an example of the spatial distribntion of electric fields and Poynting finxes in NOAA 
11158, we selected two instances before Tpre = 01 : 36 UT and after Tpost = 02 : 12 UT the 
X2.2 flare. These instances are of particnlar interest since the photospheric vector magnetic 


field changed abrnptly dnring the flare (Petrie 2012; Wang et ah 2014), and varions flare 


signatnres, snch as flare ribbons, have been observed dnring this time frame. Here we 
investigate the changes in velocities, electric fields, and Poynting finxes associated with 
these changes. We note that the 12-minnte-cadence data shown here are derived from a 
tapered temporal average that is performed every 720 seconds nsing observations collected 
over 22.5 minntes (1350 s) to rednce noise and minimize the effects of solar oscillations 
(Hoeksema et al(]2014). 


Magnetic Field: Bpreflare Bpostflare 

Fignre|^ shows the vector magnetic fields B centered at NOAA 11158 before (prefiare, 
Tpre = 01 : 36 UT) and the first time step after the onset of the flare (postfiare. 















17 


Tpost = 02 : 12 UT) and also the difference image between the two (right panel). As 
seen from the difference image on the right, the horizontal magnetic held close to PIL 
increased by over 300 G during the hare (see arrows), while the vertical magnetic held 
remained nearly constant. This horizontal-held increase is consistent with the magnetic 


held contraction scenario of Hudson et ah (2008) and is described in detail by Petrie 


(2012); Wang et ah (2014) and Sun et ah| (2012). On the diherence image we also notice 
the two circular patterns, directed counter-clockwise in negative (pixel coordinates [20,90]) 
and clockwise in positive (pixel coordinates [120,80]) polarities, meaning that the held 
connecting positive and negative polarities becomes more left-handed. Since both prehare 
and posthare magnetic helds have a right-handed twist (i.e. horizontal magnetic helds 
have clockwise orientation in the negative polarity and counter-clockwise orientation in 
the positive polarity), the observed change, dB, decreases the twist present in the prehare 
magnetogram - direct evidence of an abrupt magnetic twist decrease in both sunspots 


during the hare, which contradicts arguments (Melrose 1995) that the vertical current 
density through the photosphere should not change on the hare time scale. The sudden 
change in twist might arise from removal of magnetic helicity from the active-region’s 


magnetic held by a coronal mass ejection (Longcope and Welsch 2000 Petrie 2012). 


y^.1.2. Velocity Field. A^pi.0flare OjTld A^postflare 

For completeness, in Figure we show the three components of the velocity held 
at the times before and after the hare as well as the diherence image between the two. 
The Doppler velocity has a strong (1 km s“^) uphow close to PIL and in the sunspot’s 
umbra. During the hare the Doppler speed does not exhibit any prominent changes (see the 
background of the right panel), although such changes have been reported in other hares 
(Deng et ah 


2006). However the horizontal velocity held does change: there is an apparent 


























18 



0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140 

X-distance unit [360 km] X-distonce unit [360 km] 


Fig. 4.— Horizontal (arrows) and vertical (grayscale backgronnd) components of the mag¬ 
netic field in NOAA 11158 at preflare {top left) and postflare times {bottom left), and the 
difference horizontal held between the two, posthare minus prehare {right panel). The back¬ 
ground white and black colors represent positive and negative vertical magnetic helds, re¬ 
spectively. The blue and red colors correspond to horizontal held in areas of positive and 
negative values of the background vertical magnetic held B^. The white and black contours 
outline the positive and negative vertical magnetic helds at Bz = ±400 Gauss. The arrows 
in the right bottom corners show scales for horizontal magnetic held vectors, B/^ (left panels) 
or its change, dB/j (right panel). 


0.2 km/s drift of magnetic held away from the PIL during the hare. More details on 


the horizontal velocity pattern during the hare can be found in Wang et ah (2014), where 
45-second cadence HMI intensity maps were used. 
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X-distonce unit [360 km] X-distonce unit [360 km] 

Fig. 5.— Horizontal (arrows) and vertical (background) velocity field in NOAA 11158 at 
prefiare {top left) and postfiare times {bottom left), and the difference image between the 
two {right panel). The background white and black colors represent positive and negative 
Doppler velocities, respectively, where positive is toward the viewer (opposite to the usual 
astrophysical convention). The blue and red colors correspond to V/j originating in areas of 
positive and negative values of background vertical velocity. The white and black contours 
outline the positive and negative vertical magnetic fields at = ±400 Gauss. The arrows 
in the right bottom corners show scales for horizontal velocities. 

f.1.3. Electric Field. Epj. 0 fjaj.e ond EpQg^fj^re 

Finally, in Figure we show the vector electric field maps (“vector electrograms”) 
calculated using the PDFI method before and after the flare, and the difference between 
the two. The horizontal electric fields E/i range from 0 to 1.2 V/cm, with the strongest 


































20 


values concentrated close to the PIL. The strong E/^ at the PIL is a consequence of the 
steady velocity upflows and large horizontal magnetic helds along the PIL that lead to a 
non-inductive, horizontal electric field, — the “D” in PDFI — oriented in a perpendicular 
direction to the PIL. Given the flare-driven increase in the horizontal magnetic field at 
and near the flaring PIL, this E^ and hence E/j also increase, by up to 0.5 V/cm close 
to PIL, as shown in Figure]^ As for the vertical component, varies within [—2.2,1.8] 
V/cm. The black island at the PIL before and after the flare reflects a presence of the 
strong vertical electric held of ~ —1.5 V/cm directed into the Sun that remains nearly 
constant during the hare. What are the physical consequences of the presence of a strong 
vertical electric held? The spatial distribution of the vertical component of the inductive 
electric held from PTD, E^, is related to the time derivative of vertical current density via 
dtJz = c(V^cE^)/(47r). Hence, the presence of a nonzero E^ is related to changes in the 


vertical current, which is mostly concentrated close to the main PIL (Petrie 2013 Janvier 


et al. 2014). 


4-1.4- Poynting Flux Vector: Spreflare and Spostflare 

Figure shows snapshots of the three components of the Poynting hux vector, 

S = ^E X B, before and after the hare and also the diherence image between the two. 
The vertical Poynting varies within [—0.6, 2.3] x 10^° ergs ■ cm~‘^ ■ s~^. The maximum 
value of this range is around three times smaller than the steady-state photospheric energy 
hux estimated from the Stefan-Boltzmann law, given the temperature Tphotosphere = 5500iF, 
/ ~ 6 X 10^° erg ■ cm“^ • s“^. By analyzing the distribution of the vertical Poynting hux, 
we hnd that more than half of the spatially integrated signed Poynting hux both at prehare 
and posthare instants is injected in the range of Sz = [10®, 10^®] ergs ■cm~‘^ ■ while 
the rest of the energy is injected in the range of Sz = [10®, 10®] ergs ■cm~‘^ ■ with the 
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Fig. 6.— “Electrogram”: Horizontal (arrows) and vertical (backgronnd) electric field compo¬ 
nents in NOAA 11158 at prefiare and postfiare times, and the difference image between the 
two. The bine and red colors correspond to E/j in areas of positive and negative backgronnd 
Ez- The white and black contonrs ontline the positive and negative vertical magnetic fields 
at Bz = ±400 Ganss. The arrows in the right bottom corners show scales of horizontal 
vectors. 

positive finx. We find the strongest vertical Poynting fluxes close to the main PIL; during 
the flare the mean values of positive and negative Sz contributions within a white dashed 
box, shown on the right panel of Figure change from (1.8 to 4.5) x 10® ergs • cm“^-s“^ 
and from (—0.6 to —1.4) xlO® ergs ■ cm“^-s“^ respectively, or spatially integrating over the 
same box in terms of energy injection rate, from 1.3 x 10^^ ergs s“^ to 4.9 x 10^^ ergs s“^ 
(see white enhancement on the right panel). Summed over the 22 minutes of the GOES 
flare time for this event, the energy crossing the photosphere during the flare within the 
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box, would amount to about 4.1 x 10^° erg, that is very small compared to typical energies 
released in large flares such as this (10^^ erg or more). We also note that the absolute value 
of the spatially integrated positive Poynting flux (over NOAA 11158 ) is more that two 
times larger than the absolute value of the spatially integrated negative Poynting flux - 
more magnetic energy is transported through the photosphere upward into the corona than 
downward. 

Similar to the evolution in Eh, the magnitude of the horizontal component of the 
Poynting flux, Sh, increases near the PIL by up to 8 x 10® ergs ■ cm“^-s“^. As noted above, 
the combination of steady upflows found in the Dopplergrams at and near the PIL with the 
systematic increase in the transverse magnetic fields along the PIL leads to an increase in 
the non-inductive, horizontal electric held, E^, which in turn produces an increase in the 
Poynting flux there. 

Summarizing some of the changes in the photosphere during the flare, in Figure]^ we 
show the vertical and horizontal magnetic fields, horizontal electric held and the vertical 
Poynting flux before and after the flare and also the scans along a fixed x value (or column 
in the 2D image) through the center of the active region across the PIL at those two 
moments (right column). From the right panel we And that while the vertical magnetic held 
does not exhibit any significant changes (see i^^-panel, first row), the horizontal magnetic 
held increases at the PIL by over 300 Gauss (see i?/i-panel, second row). This leads to an 
increase in the horizontal Doppler electric held by up to 0.5 V/cm close to PIL and by 
almost 1 V/cm in some locations away from the PIL (see blue and red plots for Eh close to 
vertical dotted lines, third row) and an increase in the vertical Poynting flux close to PIL 
by up to 10^® ergs-cm“^s“^, from [1.25 to 2.25] xlO^® ergs-cm“^s“^. 
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S: 01:36 UT (preflare) dS: Postflare-Preflare 
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Fig. 7.— Horizontal (arrows) and vertical (backgronnd) Poynting vector field components 
in NOAA 11158 at prefiare and postfiare times, and the difference image between the two. 
The bine and red colors correspond to Sh in positive and negative areas of the background 
Sz- The white and black contours outline the positive and negative vertical magnetic fields 
at Bz = ±400 Gauss. The arrows in the right bottom corners show scales of horizontal 
vectors. The range of the background Sz is [—1,1] x 10^*^ ergs cm“^ s“^ in the left panels 
and [—0.4, 0.4] x 10^° erg cm“^ s“^ in the right panel. The white dashed box in the center 
of the right panel indicates the field of view where energy injection rates are calculated (see 


Section 4.1.4). 


4.2. Six-Day Evolution of Vertical Energy and Helicity Fluxes 


To quantify the long-term evolution of the energy and helicity fluxes in NOAA 11158, 
we integrate the photospheric Poynting and helicity flux maps over the AR’s field of view 
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(665 X 645-pixels) and analyze their behavior over 6 days (768 time steps). 


4 . 2 . 1 . Evolution of Free, Potential, and Total Energy Fluxes 

Panels (A-F) of Figure]^ show the six-day time evolution of area-integrated potential 
and free vertical energy fluxes and their sum, the vertical Poynting flux: [Sz^p, Szj, The 
potential component, Sz^p is calculated by taking a time derivative of the coronal potential 
energy, computed as a surface integral over the photosphere, with the potential function 
computed using a Green’s function technique. The free component is the difference between 
Sz and Sz^p. Szj = Sz — Sz^p. In the three insets on the left panels we also show evolution of 
Sz from one hour before to four hours after the flare. In panels A and C, large fluctuations 
are present in Sz,p and Sz (which were calculated independently), with timescales ranging 
from about 4 — 24 hr. In panel B, Szj does not exhibit such strong fluctuations at shorter 
timescales, which evidently cancel out in the difference Sz - Sz^p] a residual 24-hr. periodicity 
can still be seen, however. These fluctuations are related to orbital motions of the SDO 


satellite (Liu et ah 2012b) and are present over the whole held of view. Unfortunately, these 


huctuations are internal to the HMI polarization measurements, hence their removal is a 
time-consuming and complicated task. Efforts to address these artifacts are underway, but 
they have not been fully addressed at the time this article was written. Fortunately, when 
these huctuations are integrated over several hours, they do not cause large huctuations in 
the total cumulative huxes (see panel E). 

How do Sz^p, Szj, and Sz change during the hare? Insets of Panels A-C of Figure 
show that Szj and Sz increase during the hare, after which they return to nearly prehare 
values. We can understand this change in the following way. Since the change in the 


magnetic held close to the PIL is practically a step function (Petrie 2012), the value of 


the Poynting hux, which involves the temporal derivative of magnetic held, should indeed 
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appear like a delta function in time. We conclude that while the transients in the vertical 
Poynting flux include spurious signals due to flare-induced effects on the HMI spectral line 


(Maurya et ah 2012), the observed peak in the Poynting flux around the flare time is real. 
(The spatial distributions of the Poynting flux before and after the flare are shown in the 
bottom row of Figure |^) 


Panel D shows a steady growth in the region’s total unsigned flux after February 14, 
after the bulk of the region’s flux has emerged. The growth at the time of the flare is 
consistent with the steady upflows seen along the main PIT, which we expect to carry 
magnetic fields upward from the solar interior, across the photosphere. A slight 24-hr 
periodicity is discernible from 14 February onward. 


As shown in panel E, if we integrate and in time, the fluctuations 

seen in panels A-C largely disappear. Notably, no contribution from the flare transient 
is obvious. We also observe that, for a time, before February 14 12:00 UT, the active 
region appears to possess negative free energy, a spurious result. A probable explanation 
for this is that our Poynting flux estimates are affected by the orbital motions, while the 
potential energy that depends only on suffers differently from this systematic error. An 
alternative approach that might ameliorate the negative free energy budget is modification 
of boundary conditions used to compute the potential held - if instead of the Neumann 
boundary condition, given by we had used a “hybrid” Dirichlet-Neumann boundary 


condition (Welsch and Fisher 2015), which uses both (V/i ■ Bh) and B^ from the HMI data, 
our estimate of the potential energy would be Ep^hybrid = 4.2 x 10^^ erg by the hare time (cf. 
Ep = 8.6 X 10^^ erg ) implying three times more free energy Ef^hybrid = 6.4 x 10^^ erg. 


To conclude, taking into account the 29% uncertainty in Sz, due to uncertainties in 


the HMI data and the PDFI method (see Section]^), we hnd the total energy that entered 
the photosphere by the hare time to be E = [10.6 ±3.1] x 10^^ erg, which consists of 
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Ep = 8.6 X 10'^^ erg of potential energy and Ej = 2.0 x erg of free energy. 


4-2.2. Helicity Flux Evolution 


Figure 10 shows time evolution of area integrated helicity flux rates (left) and 


helicity fluxes Hr (right) in NOAA 11158 calculated using the PDFI electric helds (PDFI) 
and the velocity held alone (DFI). The apparent large-scale 24-hour periodicity is most 


likely due to signal contamination from the satellite motion (Liu et ah 2012b). Until 12:00 


UT on February 14 2011, the two helicity hux rates match each other, but at later times 
they diverge: from the velocity held (DFI, black) becomes signihcantly larger than that 

from the PDFI electric held (green). What could be the cause of this diherence? Note that 
the primary extra input to the PDFI results compared to the DFI results is information 
about Jz, which sets the inductive E^; this, in turn, is coupled to E/j through the ideal 
Ohm’s constraint, E ■ B = 0. To investigate the role that the ideal Ohm’s constraint plays 
in the PDFI helicity estimate, we calculated the helicity hux rate from the non-ideal PDF (i.e. 


PTD-Doppler-FLCT) electric helds (see blue curve on Figure 10). We hnd that the PDF 
and the DFI helicity hux rates have a similar time evolution, although the PDF helicity hux 
rate is somewhat smaller than the DFI rate. The discrepancy between PDFI and DFI helicity 
hux rates could arise because the DFI helicity hux rate is only sensitive to emergence and 
shearing. Shearing motions however, can then induce the photospheric held to unwind via 
the propagation of torsional Alfven waves into the interior, along with the accompanying 
Jz and Ez signatures, that only PDFI captures (cf., the transient twist variations modeled 


by Longcope and Welsch 2000). The constraint E ■ B = 0 then couples the resulting change 


in Ez (and Jz) to a change in Eh, thereby further reducing the PDFI helicity hux rate. 
Integrating over time we hnd that the total accumulated helicity hux injected through the 
photosphere from the start until the hare time (right panel) is Hrrpj = 7 x 10^^ Mx^, 
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which is similar to the value from PDFI; Hrp£,pi = 8.5 x 10^^ Mx^. 
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Fig. 8.— Left and Middle columns: Vertical and horizontal magnetic fields, horizontal electric 
field, and the vertical Poynting finx before and after the flare. Right column: Vertical scans 
for 5^, Bfi, Eh and Sz before (bine) and after (red) the flare at the locations shown with 
bine and red lines in the two left colnmns. Horizontal dotted lines on left and middle panels 
correspond to vertical dotted lines on the right panel. 
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Fig. 9.— Evolution of magnetic and Poynting fluxes in NOAA 11158 during six days: (A-C) 
Area-integrated potential (A), free (B) and total (C) Poynting fluxes. Corner insets show 
evolution of the same quantities one hour before and 2 hours after the flare; the X-axis is in 
hours; (D) Area integrated unsigned magnetic flux, (E) Area and time integrated free and 
potential components of the vertical Poynting flux and their sum, (F) GOES 5-minutes soft 
X-ray light curve (1-8 A channel). The vertical dotted line indicates the GOES flare peak 
time at 01:56 UT. All the quantities were calculated within the FOV shown in A-D Panels 
in Figure 
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d H r/ d t 



Hr 



Fig. 10.— Evolution of Felicity flux rates and Felicity fluxes in NOAA 11158 during six 
days: Left: Spatially integrated Felicity flux rates calculated from velocity field (DFI, 
black) and PDFI (green) and PDF (blue) electric fields; Right Total Felicity fluxes Hji, i.e. 
quantities on tFe left panel integrated in time. 
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5. DISCUSSION 


Here, we discuss how electric fields, energy fiuxes (Poynting fiuxes), and helicity fiuxes 
estimated with the PDFI method compare to earlier results. 


Using the PDFI method, we find photospheric electric field components in NOAA 11158 
whose amplitude varies from —2 V/cm to 2 V/cm. While strong horizontal electric fields 
are present in the whole active region, with the amplitude varying from 0 to 1.2 V/cm, the 
largest concentrations of strong vertical electric field are located mostly near the PIL and 
the sunspot penumbrae with the maximum amplitudes reaching 1.5 V/cm directed inward 
to the sun. As a result of magnetic field changes during the flare, the horizontal electric 
field (from the steady upward Doppler flow crossed with the horizontal magnetic field along 
the PIL) increases by up to 0.5 V/cm perpendicular to the PIL, while the vertical electric 
field at PIL remains nearly constant Ez ~ —1.5 V/cm. 


How do our estimates of the photospheric electric fields compare to results from 
earlier work? To our knowledge, there have been very few attempts to estimate all three 


components of the electric field in the photosphere. For example, using LCT, Liu and 


Zhang (2006) examined only the vertical component of inductive electric field, finding the 
maximum of Ez to be Ez ~ [0.1 — 0.2] V/cm, i.e. significantly smaller than our estimates. 
Other attempts have focused on estimating the electric field in the corona inside the 


reconnecting current sheet (RCS) (Poletto and Kopp 1986; Wang et ah 2003, 2004 Qiu 


et al. 2002, 2004 Jing et ah 2005[). For example, using a relationship between the electric 


field along the current sheet and the observable velocity and magnetic field (Priest and 


Forbes 2000), Poletto and Kop^ (1986) derived the maximum value of electric field in the 


RCS to be 2 V/cm. In a similar way, Wang et ah (2003) found the coronal electric field 


during the two-ribbon flare occurring in two stages: the coronal electric field remained near 
1 V/cm averaged over 20 minutes during the first stage, and was followed by values of 0.1 






















































32 


V/cm over the next 2 hours. Finally, in one of the latest studies on the subject, Jing et al. 


(2005) analyzed a sample of 13 two-ribbon flares, and found the maximum electric held 
inside the RCS to be in the range of 0.2 — 5 V/cm. To summarize the above, most of the 
estimates of electric helds inside the reconnecting current sheet fall within the range of 
photospheric electric helds that we hnd here for NOAA 11158. 


Using the electric and magnetic helds, we estimate the three components of the 
photospheric Poynting hux vector. We hnd the amplitude of the maximum Poynting huxes 
in the active region to be 5*2 ~ 2 x 10^° ergs-cm“^ s“^, i.e. several times smaller than the 
steady state solar luminosity, 6 x 10^° ergs-cm“^ s“^, but of the same order of magnitude. 
If we look at the values of Sz across the AR, where most of the magnetic energy is coming 
from, we hnd even smaller values, ranging from 10® to 10^® ergs-cm“® s“^. A question 
naturally follows this calculation: Is such a Poynting hux consistent with the amount of 
magnetic energy stored in the coronal part of the active region? This is the question we 
address in Table [H 


In Table [T| we compare estimates for coronal energy and helicity for NOAA 11158 that 
we derive in this paper with results from other papers. In the hrst four rows, we show 


results from Evolutionary estimates, i.e. works by Tarr et al. (2013); Liu and Schuck (2012); 


Tziotziou et al. (2013); Vemareddy et al. (2012b) and Jing et al. (2012), where energy or 


helicity or both are found from the evolution of photospheric magnetic helds, by summing 
the energy or helicity hux rates injected through the photosphere from the beginning of the 
magnetogram sequence until the moment before the hare. In the last hve rows, we show 


results from Instantaneous estimates, i.e. works by Malanushenko et al. (2014); Aschwanden 


et al. (2014); Sun et al. (2012); Tziotziou et al. (2013) and Jing et al. (2012), where a single 


close to the hare time magnetogram or EUV image is used to estimate the energy or helicity 
of the corona. The “Method” column shows the type of method used to hnd the energy 
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and helicity estimates on the right. In some papers, other methods were nsed to calcnlate 
the potential energy; they are indicated with a letter next to the estimate. In addition, we 
indicate the type of inpnt data nsed in the calcnlations. This helps to explain the difference 
in resnlts between some papers that nsed the same methods, bnt applied them to different 


inpnt data, e.g. differences in estimates of the potential field energy, Ep, by Tarr et ah 


(2013) and Snn et al. (2012), where Blos and have been nsed respectively. 


The total magnetic energy from different models, shown in Table ranges from 
6 X 10^^ ergs (Malannshenko et al. 2014) to 12 x 10^^ ergs (Tziotzion et al. 2013). We 


noticed that both coronal NLFFF methods, which nse the EUV coronal loops instead of the 


transverse magnetic field as a constraint for the NLFFF extrapolation (Malannshenko et al. 


2014 Aschwanden et al. 2014), derive total energies that are the smallest of all total energy 
estimates: 6 x 10^^ ergs and 8.6 x 10^^ ergs respectively. In contrast, the photospheric 


NLFFF methods, which nse a vector magnetogram for extrapolation (Snn et al. 2012 


Tziotzion et al. 20131), yield the largest estimates for the total energy: 10.6 x 10^^ ergs 


and 12 x 10^^ ergs respectively. The evolntionary estimates, that derive the total energy 
by integrating the energy flnx, inferred from the photospheric velocity or electric fields. 


yield total energies in between the two. For example, nsing the MCC method, Tarr et ah 
(2013) derive E = 8.5 x lO^^ergs. Using DAVE4VM approach, Lin et al.l (2012b) and 


Tziotzion et al.l (120131) yield similar estimates oi E = 8.8 x 10^^ ergs and E = 8.0 x 10 
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ergs, respectively. Finally, in this paper nsing the PDFI method we hnd the total energy 
of E = [10.6 ± 3.1] X 10^^ ergs. To conclnde, taking the PDFI and HMI nncertainties into 
acconnt, we find that the total energy, estimated right before the flare, is consistent with 
E from DAVE4VM, MCC and NLFFF and is slightly larger than the coronal NLFFF 
estimates. 


If we look at the temporal evolntion of energy A, we notice that E derived from PDFI 
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(Figure is consistent with E from DAVE4VM (Figure 14 in Liu et al. ( 2012b[ )) and 


NLFFF (Figure 4 in Sun et al. (2012)). In fact, the energies E derived from the PDFI and 


DAVE4VM are almost identical until February 14 18:00, and then start diverging several 
hours before the flare. By the end of magnetogram sequence, at 18:00 UT on 16 February, 
the E from DAVE4VM and PDFI are 12 x 10^^ ergs and 14 x 10^^ ergs, respectively. 
Still, this discrepancy lies within the uncertainty of the total energy estimate (29%, see 


Section 3.3). Comparing the temporal evolution of E from PDFI and NLFFFs, we again 


End that they are consistent with each other, with less than 15% differences between the 


two, which is within our uncertainty (29%, see Section 3.3). 


The potential field energy from different models, also shown in Table ranges from 


4.8 X 10^^ ergs (Malanushenko et al. 2014) to 8.6 x 10^^ ergs (this paper). Similar to our 


approach. Sun et al. (2012) used the Green’s function and to estimate potential held 


energy of 8.0 x 10^^ ergs, that, within the HMI uncertainty, is consistent with our estimate. 


The potential held energies calculated by Tarr et al. (2013); Malanushenko et al. (2014) 


and Aschwanden et al. (2014) are smaller than Ep from this paper and from Sun et al. 


(2012). This diherence might be due to the fact, that in contrast to this paper and Sun 


et al. 

(2012 

), where B^ is used to calculate Ep, 

Tarr et al. 

(2013 

); 

Malanushenko et al. 


(2014) and Aschwanden et al. (2014) used B^os that tend to underestimate the held at 


Blos > 1000 G due to limitation of the LOS pipeline algorithm (see Section 7.2.3 and 


Figure 17 in Hoeksema et al. 2014) and therefore lead to underestimate of energies by a 


factor of 1.4 — 1.7 ( Liu et ah]|2012bt [Malanushenko et al.||2014 ), that we hnd to be consistent 
with the diherences in Ep in Table 

The free magnetic energy from diherent models, also shown in Table ranges from 


1.0 X 10'^^ ergs (Aschwanden et al. 2014) to 2.9 x lO'^^ ergs (Tarr et al. 2013). Using the 


MCC model and the hare ribbon locations, which allow one to derive the footprint of the 
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reconnecting magnetic fields, Tarr et al. (2013) find the initial pre-flare free energy in the 


corona to be Ef = 2.9 x 10^^ ergs, consistent with estimates of Sun et al. (2012) and 


the results of this paper. Tarr et al. (2013) also hnd that more than 50% of this energy, 
dE = 1.7 X 10^^ ergs, is released during the flare. Coronal NLFFF methods (Aschwanden 


et al. 2014; Malanushenko et al. 2014) hnd that 60% to 80% of the free energy is released 


during the hare, but the values of the pre-hare free energy Ef that they hnd are roughly 


factors of two (or more) times smaller than the values from Tarr et al. (2013); Sun et al 


(2012), and this paper. 


Another important fact one must keep in mind when comparing cumulative free 
Poynting huxes and coronal free energies, is that the total Poynting hux only gives us 
information about the total energy that entered the corona from the photosphere. What 
we do not know is how much of this total energy leaves the corona from above into the 
heliosphere during eruptive hares prior the X2.2 hare. For this reason, our PDFI free 
energy estimate is the upper limit of the energy available in the corona. To summarize, the 
free magnetic energy, which we hnd from a diherence of integrated Poynting hux and the 
potential energy, is up to two times larger than the free energy estimated from the coronal 
NLFFF codes, and 20 — 30% smaller than free-energy estimates from the MCC model 
and photospheric NLFFF extrapolations (Tarr et al. [2013 ; Sun et al. 2012). The temporal 
evolution of the PDFI free energy is within 20 — 30% of Ef{t) from the photospheric NLFFF 
code. 


Finally, in the last column of Table [T| we compare the total relative magnetic helicities 
before the X2.2 hare. The time integrated total helicity hux we hnd with DFI and PDFI 
techniques is = 7 x 10‘^^Ma;^ and Hppppj = 8.5 x 10^‘^Mx'^ respectively. The 

diherence between the two is consistent with the results we found for the ANMHD test 
case, where the DFI method reconstructed the total helicity around 10% more accurately 
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than the PDFI method (see Table 3 in Kazachenko et ah (2014)). For comparison, nsing 


DAVE4VM, [Lin et ah (2012b) and Tziotzion et ah (2013) fonnd that the total amonnt 
of helicity injected into corona is 6.5 x 10^^ Mx^ and 8.5 x 10^^ Mx^ respectively. These 
resnlts are consistent with onr PDFI estimate, given the differences between the DAVE4VM 


and the PDFI accnracies (Kazachenko et al. 2014) - the ratio between the total PDFI 


and DAVE4VM-reconstrncted helicity flnxes and the actnal helicity flnx are 0.94 and 1.0 


respectively for the ANMHD test case (see Table 4 in Kazachenko et al. (2014)). Using a 


different velocity reconstrnction method, DAVE, Vemareddy et al. (2012b) and Jing et al 


(2012) find the total relative helicity of 6.0 x 10^^ Mx"^ and 5.5 x lO"^"^ Mx^ respectively. As 


shown by Schnck (2008) for the ANMHD test case, combining DAVE flows with ANMHD’s 


vertical velocity overestimates the total helicity flnx by at least 40%, hence the disagreement 
between DAVE and PDFI resnlts for NOAA 11158 is not snrprising. Finally, nsing different 
NLFFF approaches, Tziotzion et al. (2013) and Jing et ^ (2012) End Hr = 13 x 10^^ Mx^ 
and Hr = 5.2 x 10^^ Mx^ respectively. To snmmarize, nsing the PDFI method we find the 


relative magnetic helicity consistent with DAVE4VM estimates (Lin et al. 2012b Tziotzion 


et al.||2013), bnt very different from (and in between) the NLFFF estimates (Tziotzion et al. 


2013 Jing et al. 2012). 


6. CONCLUSION 

The electric held on the Snn plays an important role in transporting energy, heating 
plasma, and accelerating and transporting charged particles. Estimates of photospheric 
electric and magnetic held vectors make the estimation of Poynting flnx of electromagnetic 
energy crossing the photosphere and the flnx of relative magnetic helicity straightforward. 
Taking advantage of the newly released high temporal and spatial resolntion HMI vector 


magnetograms (Schon et al. 


2012), and the recently developed PDFI electric-field inversion 
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Table 1: Summary Table: Comparison of coronal energy and lielicity in NOAA 11158 where 
an X2.2 flare occurred on February 15 2011 01:35 UT: dE - change in the coronal free 
magnetic energy during the flare, [Ef, Ep, E] - free, potential and total magnetic energies 
in the corona before the flare at 01:35 UT, Hji - relative magnetic helicity in the corona 
before the flare at 01:35 UT. In Evolutionary or Photospheric estimates, coronal energy and 
helicity are calculated cumulatively from tracking magnetic held evolution, i.e. by integrating 
photospheric energy and helicity hux rates in time. In Instantaneous or Coronal Estimates 
energy and helicity are calculated instantaneously, using extrapolation of the photospheric 
vector magnetic helds. All quantities have been calculated using the indicated Method, unless 
a more specihc method is used, like in the case of potential energies (see letters U)-(/)) 
photospheric estimates we use a start time of 14:11 UT on 10 February. 
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method (Kazachenko et al. 2014), we apply the method to a six-day sequence of vector 
magnetic held measurements of NOAA 11158, from 10 to 16 February 2011. From these 
12-minute cadence measurements, we derive the evolution of electric helds, the Poynting 
hux, and the helicity hux. During the interval of study, an X2.2 hare occurred, along with 
35 M- and C-class hares. 


We analyze the spatial distribution of the derived electric held and Poynting hux 
maps, their temporal evolution, and changes during the X2.2 hare. We compare derived 
PDFI electric helds with various estimates of a typical coronal and photospheric electric 
helds made to date, PDFI energies and helicities with those previously reported in the 
literature for this active region (NLFFF, MCC and DAVE4VM estimates). The results are 
the following: 


1. We hnd the photospheric electric held vector, which typically ranges from —2 to 2 
V/cm, to increase its magnitude by up to 0.5 V/cm at the PIL and 1 V/cm away 
from the PIL during the hare. The horizontal component is mostly concentrated 
along the PIL, while the vertical component is largest at the PIL and in the sunspots’ 
penumbrae. The range of photospheric electric helds is consistent with the coronal 
electric helds, derived from the motion of hare ribbons. 

2. We hnd the photospheric Poynting hux ranging from [—0.6 to 2] x 10^° erg-cm“^s“^ 
with majority of the energy hux moving upward into corona and more than half of the 
total energy input rate injected from within the range of [10® to 10^®] erg-cm“®s“^. 
The largest vertical Poynting hux is concentrated at the PIL. 

3. Integrating the Poynting hux in time we hnd the total magnetic energy before the 
hare, E = [10.6 ± 3.1] x 10^^ ergs. In spite of a very diherent approach, it is consistent 
within the uncertainty with the total energies from DAVE4VM, MCC and NLFFF 
methods’ estimates and larger than the coronal NLFFF estimates. 
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4. The potential field magnetic energy before the flare, estimated via the Green’s 
fnnction, Ep = 8.6 x 10^^ ergs, is consistent with Ep = 8.0 x lO^^ergs (Snn et al. 2012) 


Tarr et al. (2013); Malanushenko et al. (2014) and Aschwanden et al. (2014) derive 


a similar or somewhat smaller Ep, using different computation methods and Blos 
instead of that is prone to being underestimated in the strong field regions due to 
limitation of the LOS pipeline algorithm (Hoeksema et al. 2014|). 


5. The free magnetic energy before the flare from the PDFI method, Ef = 2.0 x 10^^ 
ergs, is up to two times larger than the free energy estimated by coronal NLFFF 
codes ( Aschwanden et ^|2014 Malanushenko et al. 2014), and 20 — 30% smaller than 
the free energy estimates from MCC model and photospheric NLFFF extrapolations 


(Tarr et al. 2013 Sun et al. 


6. Analyzing the temporal evolution of cumulative energy, E, from PDFI, NLFFFs and 
DAVE4VM, we find less than 15% differences between the three. The temporal 
evolution of the PDFI Ef is less than 20 — 30% different from the Ej from the 
photospheric NLFFF codes and several times larger than Ef from the coronal NLFFF 


codes (Aschwanden et al. 1120 14; Malanushenko et al. 2014). 


7. We find the relative magnetic helicity to be consistent with DAVE4VM estimates 


(Liu et al. 2012b Tziotziou et al. 2013), but very different from the NLFFF estimates 


(Tziotziou et al. 2013 Jing et ar||2012 ). 


8. Using Monte-Carlo simulations, we find that the levels of the errors in the HMI data 
lead to uncertainties in the horizontal electric field of 13% to 18% that result in 
errors in the vertical Poynting flux of around 14%. If we add those uncertainties to 


the errors that we found when testing the PDFI method (< 25% in Sz, Kazachenko 


et al. (2014)), then we end up with 14% to 29% errors in the vertical Poynting flux 


depending on the LOS angle. Also, since some free energy might have been released 
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by flares prior to the X2.2 flare, we view the PDFI estimates of the free energy injected 
through the photosphere to be upper limits on the total free energy available before 
the flare. 


This study is the first application of the PDFI electric field inversion technique to 
photospheric vector magnetic field and Doppler measurements. We find that the total 
amount of energy and helicity injected through the photosphere before the flare estimated 
by the PDFI method is consistent with estimates from other approaches, in spite of differing 
techniques. This agreement is very promising, implying that the PDFI technique is not 
only capable of describing the coronal energy and helicity budget, but can also provide 
instantaneous estimates of energy and helicity transferred through the photosphere. 


We believe that both the derived dataset of PDFI electric fields and the PDFI 
method itself will be useful to the science community for analysis of the evolution and 
spatial distribution of the photospheric electric helds, fluxes of energy and helicity, and 
their relationships with flare activity. In addition, PDFI electric helds can be used as 
time-dependent boundary conditions for data-driven models of coronal magnetic held 


evolution (Fisher et ah 2015). The dataset for NOAA 11158 is available for downloading on 
our websit4El 


^http://cgem.Stanford.edu 
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A. SCALING B AND E TO CARTESIAN MERCATOR MESH 

The data we analyze have been transformed from plane-of-sky to Mercator projection 
with a local, Cartesian coordinate system centered on NOAA 11158. The distortion of pixel 
scales in this transformation has implications for inferring electric fields from magnetic 
evolntion that we describe in detail here. 

When observed magnetic fields are reprojected from the Snn’s observed plane-of-sky 
snrface onto a plane nsing a conformal (angle-preserving) mapping, pixel areas in different 
regions of the new coordinate system generally do not correspond to the same physical areas 
on the Snn. Denoting the area of a pixel ij on the Snn as Afj, and in the mapped coordinates 
as Afj , Afj ^ Afj near the center of the projection, bnt Afj ^ Afj far from center. In 
contrast, the reprojection does not directly alter the magnetic field valnes themselves: the 
fields are interpolated onto the new grid point, bnt the interpolation attempts to faithfnlly 
represent measnred held valnes at each point. While observed magnetic held vector data 
have nnits of held strength, i.e. Ganss, they are more accnrately described as measnrements 
of pixel-averaged hnx densities, Mx/cm^. This is becanse the hnx in pixel ij conld be 
conhned within a snbregion of the pixel with fraction area fijAfj, where the hll fraction fij 
obeys 0 < fij < 1, implying a true held strength larger that the apparent 

held strength (“pixel-averaged hnx density”) measnred by the instrnment. 

The qnestion, however, arises: Do we need to modify the held strengths in the new 
projection - call the original held strengths B'^, and the reprojected held strengths B^ - 
to acconnt for the distorted areas? If so, how? Do vertical and horizontal magnetic held 
components need to be compensated in the same way? How shonld velocities and electric 
helds be modihed to compensate for this distortion? 


Below we consider the transformation nsed in this paper, from the plane-of-sky to 
Mercator projection (Welsch et al.[2009). Briehy, the distortion of scale in cylindrical 
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projections such as Mercator is a function of latitude alone. The horizontal Mercator 
coordinate, x, is mapped one-to-one with the heliocentric longitude, 0, and is independent of 
latitude. Because the distance between lines of constant longitude decreases with increasing 
latitude, the projection’s scale (the distance of the sphere corresponding to a hxed distance 
in the projected image) must decrease (fewer solar Mm per Mercator pixel) with increasing 
latitude. As lines of constant longitude converge towards the poles on a spherical surface, 
the physical distance between such lines goes to zero. In Mercator coordinates, however, the 
distance between two lines of constant longitude Ax, is hxed and is independent of latitude 
6. Effectively, this means the projection magnihes distances towards the poles. (This 
effect is easily seen on Mercator projections of the Earth’s surface, on which Antarctica 
and Greenland, for example, appear too large relative to landmasses at lower latitudes.) 
Consequently, displacements in the vertical Mercator coordinate, dy, corresponding to a 
hxed latitudinal displacement d9, should increase towards the poles. Because the physical 
length corresponding to a hxed dy shrinks towards the poles as l/cos{6), and dx scales in 
the same way, pixel areas in the reprojected system scale as cos^{d) compared to areas on 
the Sun. 

Consider two reprojected pixels with the same normal magnetic held, one at a high 
latitude and one at a low latitude, denoting these respectively. These have the 

same huxes in the projection, ^ in each reprojected pixel, however, 

corresponds to diherent huxes on the Sun: ^ ^ ^m,h gjj^gg ^j^g actual 

solar area corresponding to the high-latitude pixels is smaller. We choose to handle this by 
multiplying Bn in each reprojected pixel by cos^(0), to compensate for distortion of its area 
as a function of latitude; at the same time, we do not rescale pixel lengths. 

When hux emerges into a pixel, transports horizontal held B^ along a pixel edge of 
length L over a time interval At, meaning = VnLBhAt has emerged (see discussion in 
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Section 2, and Fignre 3 of Welsch et al. (2013)). This means that the normal flnx in the 
pixels that share this edge must change by ^em, to account for the emerged flux. We have 
assumed that flux in each pixel sharing L has been compensated by a factor of cos‘^{6). For 
the emerged horizontal flux to match the changes in the vertical flux, two factors of cos{6) 
must be present in the product VnB^, since we are not rescaling L. We can scale by 
cos“(6') and Bh by cos^{9), which yields 


a+ 13 = 2. (Al) 

When flux is horizontally transported from one pixel to another, Vh transports the 
vertical held Bn across a pixel edge of length L over a time interval At, meaning a hux 
Export = VhLBnAt has been moved. This means the normal huxes in the pixels that share 
this edge must each change by Export, to account for the transported hux. Because Bn has 
already been rescaled by cos‘^9, Vh in the product VhLBn does not need any scaling. 

The pixel-integrated vertical Poynting hux = {E^ x Bj,) at high latitude must be 
scaled by cos‘^{9), to account for the area distortion. Recall that cEh is proportional to the 
sum of VnBh, scaled by cos°‘~^^{9), and VhBn, are already scaled by cos‘^{9). Because Sz must 
scale as cos‘^{9) for any v, we can consider the special case Vh = 0, implying Sz = VnB^, so 


a+ 213 = 2. 


(A2) 


Comparing equations ^ and |A2| yields /3 = 0, so Bh is not changed, but a = 2, so Vn 
is scaled by cos^{9). Physically, this implies that the vertical transport of emerging hux 
at high latitudes is scaled to make sure a “scaled amount” of hux emerges. This scaling 
implies that Eh is automatically and implicitly scaled (via scaling applied to Bh and n„) by 
cos‘^{9), and that Ez = Vh x Bh is unsealed. 
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Summarizing the above, we scale Bn,Vn,Eh by cos'^{9) and do not scale Bh^Vh^E^. 
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